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Abstract. We study the principal bifurcation curve of a third order equation 
which describes the nonlinear evolution of several systems with a long-wavelength 
instability. We show that the main bifurcation branch can be derived from a 
variational principle. This allows to obtain a close estimate of the complete 
branch. In particular, when the bifurcation is subcritical, the large amplitude 
stable branch can be found in a simple manner. 
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1. Introduction 

Long wavelength instabilities have been found in diverse physical systems, such as 
Rayleigh-Benard convection between insulating boundaries [I], Marangoni convection 
in the presence of a free surface 0E1E], directional solidification [HlEj, Langmuir 
circulation in the ocean |Jj and others [B]- Evolution equations that describe the 
nonlinear evolution of the instability can be derived by means of an expansion 
procedure based on the long scale of the instability 9 . Unlike the case of instabilities 
with non zero wavelength where the weakly nonlinear evolution of the instability is 
governed by the Ginzburg-Landau equation, for long wavelength instabilities different 
equations arise according to the problem considered [UJ- In the absence of a 
unique equation to describe these instabilities each equation needs to be examined 
to determine whether it provides an adequate description of the problem under study. 
It is often found that the bifurcation from the basic steady state is subcritical, thus 
unstable. In such cases it is important to determine whether large amplitude stable 
solutions exist. The usual procedure to capture the turning back of a subcritical branch 
is to perform an asymptotic expansion close to the bifurcation point. Suitable scaling 
of parameters is necessary in order to capture stable solutions. If this procedure fails 
one may search for approximate large amplitude solutions of the equation which is 
not always feasible. The alternative is to obtain the bifurcation curve numerically. 

The purpose of this work is to study a prototype long wave evolution equation 
which, in a given parameter range, bifurcates subcritically from the basic state |lllll2j . 
and show that the principal bifurcation branch derives from an integral variational 
principle. It is then possible to obtain lower bounds on the complete branch that 
are as close as desired to the exact solution. While we treat a specific equation, the 
method can be applied to other equations as well. In Section II the main result is 
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derived. In Section III a particular exactly solvable case is discussed. Analytical and 
numerical results for the general case are given in Section IV. 

2. Formulation of the problem 

We consider the long-wavelength equation 

/* + Afe + - + M/e/«)e + = °> W 

in a periodic box of length 2£. The boundary conditions f^(±L) = f^c(±L) = are 
chosen to satisfy physical requirements. With these boundary conditions the space 
average / is conserved. 

This equation describes the weakly nonlinear evolution of Rayleigh-Benard and 
Marangoni convection between insulating boundaries |1 II 1121 ITHl 1141 115| . Particular 
cases of steady state solutions to this equation have been studied in connection with 
directional solidification. In this context the case v = is known as the (steady 
state) Riley-Davis equation with vanishing segregation coefficient. The case v = 
and fi = is known as the Sivashinsky equation. The bifurcation structure of these 
two equations has been studied extensively [SJ El E3 ■ In the context of Rayleigh- 
Benard convection, the term proportional to fi arises due to the asymmetry between 
upper an lower mechanical boundary conditions [H] and the term proportional to 5 
arises when non-Boussinesq effects are included |lllll2lH4*l . The purpose of this work 
is to study the principal (monotonic) branch that bifurcates from the first eigenvalue. 
These solutions are sometimes called simvlefiT^. 

Introducing a new length scale x = £/L, scaling the amplitude as / = \TkJv 4>, 
and integrating once, the equation can be written as 

r(f> x + 4>xxx - 4>l + A(f)x<f>xx + B<f>4> x = 0, (2) 

for appropriate constants A and B. The constant of integration in the equation above 
vanishes due to the boundary conditions. Due to the symmetry x — > —x we may 
consider the solution on half the interval and impose 

<Mo) = M 1 ) = o- (3) 

In addition, for the principal branch we impose 

0(0) =0 m , 0(1) = 0, X <O. (4) 

Previous studies of this equation show that the quadratic term B<jxf> x is destabilizing, 
it gives rise to a subcritical instability (|3J 111! I12p. If B < the bifurcation is 
supercritical. The bending back of the unstable branch can be described close to 
the transition point B — by means of small amplitude expansions. Expanding the 
solution in a small parameter e, as 4>{x) = e((f>o + e<$>\ + ■■■) and the eigenvalue as 
r = ro + tT\ + we obtain for the eigenfunction in leading order 

, / s , 1 + cos [to] 
<Po{x) = <p m ^ ' 

and at order e 

The eigenvalue up to order e 2 is given by 

r = tt 2 - eB^ + e 2 ^-[187r 4 + (B — V2n 2 A) 2 - (3 - 2V2)n 2 AB}. 
2 96t: z 
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For positive B the bifurcation is subcritical, but the order e 2 term shows that it 
regains stability for small enough B. This result is valid only for small amplitude and 
sufficiently low values of B. 

Our purpose is to show that the principal branch r(0 m ) can be obtained from a 
variational principle from which lower bounds on the bifurcation curve are obtained. 
This enables us to obtain non perturbative results valid at arbitrary large amplitude 
and for all parameter values. The method used to derive this result is related to that 
developed to study a class of reaction diffusion problems |181 1191 12U| and others [21] • 

The first step is to rewrite the equation in phase space, therefore we define 

ax 

and write the problem to be solved as ^2] 

r P + P-JI (p-JI J - P 3 + Ap 2 ^ + B(j) P = 0, (5) 



with 

p(O)=p(0 m )=O, P >0. (6) 

Even though the equation above can be reduced to quadratures, an explicit expression 
for the eigenvalue r(<p m ) cannot be obtained except in the case B — (Appendix A). 
From the variational principle that we construct below, the exact solution for the case 
B = can be recovered, and arbitrarily accurate bounds can be obtained when B ^ 0. 

Let g((j>) be an arbitrary positive continuous function, with continuous second 
derivatives, such that g(0) = g(4> m ) — 0. Multiplying Eq.© by g/p and integrating 
by parts we obtain, after regrouping terms, 

r / *" g(cj>) d<f> + B [ m dtjy = [ " p 2 ^)h(4>) ty, (7) 



^0 JO Jo 

where we have defined 

Kcp) = -\g" + g + ^g'. (8) 

Boundary terms which appear when integrating by parts vanish because of © 
and the conditions on g. In addition we require that the function g be such that h > 
in (0, m ), a condition which can always be met (see Appendix B). 

We now use Holder's inequality, 

i/fc 

Mo) '<!■•) y\\ M>»l''-«'<>) / .s(,>), /( o),/„. 

with 



to bound the right side of equality Q. Choosing s = (p 2 /),) 1 / 3 , q = p 2 / 3 , I = 3 and 
fc = 3/2 we obtain 





h 1/3 d(f>. 



Using 



l 

dx 
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the above inequality becomes 



p 2 hd(j) > 



h 1/3 d(j) 



Inserting @ in (7J we obtain the desired lower bound, 



r > 



ff" &l> -Bf* m <f>g(<P)d<l> 



(9) 



(10) 



To show that this is a variational principle we must prove that there exists a function 
g = g for which equality holds in (10) ^Hl- In Holder's inequality, the case of equality 



corresponds to s 



(f 



h 



In the present case we have that equality in Q holds when 



1 

P 3 



From the definition of h we have equality in (|10|l when g is a solution of 
1 „ 



\9 



A , 1 

~ 9 = ^< 



(11) 



(12) 



2- r 

satisfying g(0) — g(<fi m ) = 0. To prove that the solution to this boundary value 
problem, which we call g, exists and that it yields finite values for the appropriate 
integrals in Ijl0|l is straightforward (see Appendix B for details). 

Since equality holds for a special function g we conclude that the principal 
bifurcation curve of is given by 

3 



max 

9 



' (13) 
where the maximum is taken among all positive continuous functions g (with 
continuous second derivatives and such that g(0) — g(4> m ) = 0) for which h > and 
integrals exist. The maximum is achieved for g = g up to a multiplicative constant. 

It is worth pointing out that the Euler-Lagrange equation for (|13fl can be obtained 
easily (Appendix C) and is, effectively, equation JSJ when g = g. 

Expressions IjlOfl and (|13fl constitute our main result. Using adequate trial 
functions g in (|10f> we obtain lower bounds on the bifurcation branch r(<j) m ). Moreover 
we know there exists a trial function for which equality can be attained. 

In the following section we study the case B — and show that the maximizing 
g can be constructed explicitly in this case. 



3. The case B=0 

In this section we show that the exact formula for the eigenvalue can be obtained from 
the variational principle when B = 0. Furthermore, the optimal g obtained for this 
case will provide a good trial function for B ^ 0. 

To obtain the solution for the eigenvalue, instead of solving the maximizing 
problem (|13|) . we address the equivalent problem of minimizing its inverse. Let 

r<l>n 
JO 



min J[g] 

g 



mm 

9 



g{4>) d(j) 



(14) 
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with g > 0, and ft, > as denned in (JHJ. Given h, g can be found in terms of h by 
means of the Green's function G(<j>, cj>') for ijSJ). Defining 

= h 1 ^^) (15) 

we have that 

g(4>) = [ m G(<j>,<t>')s 3 (0') dtf, (16) 
Jo 

where the Green's function G(0, 0') corresponding to equation JSJl is given by 

rU jj\ = J bs iJ(H m f A[ ^' )l2 Sinh (^) - 0')) for < ^ 

1 fc sinh 2 Wm ) eA(0 - 0,)/2 sinh W) sinh(6(0 m - <P)) for^>0', 



6= ^^2+8. 



with 



Then, 



where 

F{ct>)= I m G(C,^ = l-e-^/ 2 
Jo 

x cosh(60) - coth(60 m ) sinh(60) + e j40m/2 cosech(60 m ) sinh(60) . (19) 
Finally, (|14|l can be written as, 

J TO = mm J[g\ = mm — ^— , (20) 



which is a simple variational problem for which it is straightforward to find the 
minimizing function. The minimizing s is given by 

s - w = ^vwy (21) 

where 

N^f^^L. (22) 

Jo vnf) 

Using this function in l|20() . we find 

■^min ^( s min) ~j^2 

and the exact value of r for B ~ is 1/J(s m ; n ). Spelled out explicitly, 

which is the exact expression for the bifurcating branch when B=0 (Appendix A). 
Notice that we have also obtained implicitly the solution to the equation itself, since 
equality in (|T3|) holds when h = s^ in = 1/p 3 . The solution of the equation is obtained 
by integrating p = —defy/dx = l/s m j n . 
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4. Bounds for B ^ 

The exact solution for the principal branch cannot be obtained in closed form when 
B^O. Instead, upper bounds can be obtained using different trial functions in l(TU|) . 
In this section accurate bounds will be obtained numerically. A simple analytical 
bound will be obtained in the special case A = 0. 



4--1- Numerical bounds 

Guided by the results of the previous section, we choose as a trial function the exact 
solution for the case B=0, that is, 

gxffl = T m G(0,0O4 in ( (23) 



Replacing this choice for g in (|ll)fl we obtain 

2 



r > 




/n 



)-fr^(r* w '*' w )* (24) 



s/FW). 
where N is defined in 122jl . 

A different, simpler trial function, is obtained choosing h(4>)= 1, which gives 

52 (0) = [ m G(0,0')#'- (25) 



n 



In Figure 1 below we show the bifurcation curve for A=0, B=5. The solid line is 
the exact numerical solution obtained using the software Auto (22- The dot-dashed 
line is the lower bound obtained using the trial function 1 2. 'ill , and the dashed line is the 
bound obtained with the trial function l|25[) . The bound obtained using gi as a trial 
function provides a very close estimation of the bifurcating branch, especially at low 
amplitudes. Already the simple trial function 52 which does not give a close estimate 
is enough to prove that the subcritical bifurcation is stabilized at larger amplitudes. 
In Figure 2 we show the numerical results for A=0, B=10. The curves shown are 
as in Figure 1. At small amplitudes g\ gives a close estimate, at large amplitudes it 
captures the qualitative behavior, still providing a good bound. Here too the simple 
trial function 52 suffices to capture the stabilization of the subcritical branch at larger 
amplitudes. 

In Figure 3 we show the numerical results for A=10, B=10. The solid line is the 
numerical solution of the bifurcation problem using Auto, the dashed line is the bound 
obtained with g\. 

4-2. Analytical bounds 

In the case A=0 a simple admissible trial function is g^{(t>) — (t>{4> m — </>)■ Then, 

h(<j>) = 1 + (j)(<f> m ~4>)> <t>{<j>m - <t>)- 
Using gy,{4>) and this lower bound for the corresponding h((f>) in (|10|) we obtain 

^ [4>{(p m - 0)]V3 d<p — B f* m 2 (0 m - 0) # 



r > 



fa m <H0 m - 4>) • 



> 65(4/3, 4/3)0, 2 „-|0 m 

>3,18^-|0 ro , (26) 
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Figure 1. Exact bifurcation curve and bounds with different trial functions. 




rln 2 

Figure 2. As in Fig.l for different parameter values. 



where B(x, y) is the standard Beta function. This bound, valid for all amplitudes and 
all values of B shows that the subcritical branch always turns around and becomes 
stable at sufficiently large amplitude. 

4-3. Bounds for solutions of zero average 

Here we address the problem of solutions to of zero average. The bifurcating branch 
f{4>m) obtained above does not correspond to a solution of fixed average. Each point 
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Figure 3. Exact bifurcation curve and bound with the B=0 exact solution 



along the branch has different average, <\> = 4>(4> m )- It is possible to obtain bounds for 
solutions of zero average as follows. Define u — 4> — cj>. In terms of u, equation 
becomes 

(r + B<j>)u x + u xxx - u x + Au x u xx + Buu x = 0, 

with 

«a:(0) = U X {1) = 0, u(0) = <p m - </>, u{\) = -<t>. 

The amplitude of the function u is <fi m , its average is zero and the bifurcating branch 
ru =0 is given by 

r%=o(0 TO ) = r(<p m ) + B(f>. 

Since < <f> < <fi m , we know that for B > 0, B<fi > 0, whereas for B < 0, Bcj) > B<f> m . 
Then we know that for a solution of zero average 

> / r(<f> m ) for B > , . 

r " =0 " \ r(4> m ) + B<\> m for B < 0. 

The analytical bound (|26|) for A — becomes for solutions of zero average, 

f 3.18^ - S0 m /2 for B > 
ril=0 ^ \ 3.18^ + B^„/2 for B < 0. 

This shows that also for zero average solutions, all subcritical branches turn back at 
sufficiently large amplitude. 

5. Summary 

We have studied the principal bifurcation branch of a third order equation which 
exhibits a subcritical bifurcation for certain parameter values. The bifurcation 
problem can be characterized by an integral variational principle from which lower 
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bounds on the complete branch can be found. With the use of appropriate trial 
functions it is possible to obtain an accurate estimation of the complete branch, 
without resource to perturbation theory. This permits to capture the stabilization 
of the subcritical branch even when the turning point occurs at large amplitude. 

While we focused on a specific equation, the method can be applied to other 
nonlinear equations [18-21]. 
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Appendix A. 

In this appendix we show that even though the problem the problem can be reduced 
to quadratures, an explicit expression for the eigenvalue can be obtained only when 



The solution to this equation which satisfies the boundary conditions it(0) = u(<f> m ) — 
is given by 



B = 0. 



Equation (J^J can be reduced to a linear equation for u = p 2 , 




(A.l) 



u = {r + —)F{^)+BZ{(t>), 



(A.2) 




cosech(60 m ) sinh(60). 



(A.3) 
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The eigenvalue r follows from the condition 

' * d<t> 



lo p{4>) 

which may be written as 



1, (A.4) 



= 1. (A.5) 



lo ^{r + AB/2)F{<P)+BZ{<P) 

Since this integral cannot be performed analytically, only in the case B — an explicit 
expression for r is found. 

Appendix B. 

In this appendix we analyse the properties of the solution of <|12fl . i.e., the solution of 
the two point boundary value problem 

- 1 - 9 » + ± 9 > +9 =± (B.l) 

with g(0) = g((f> m ) — 0, where p(<f>) is a solution of (J5J with p(0) — p(4> m ) — 0. One 
can conveniently write l|B.l[) as 

L{g) = ^e- A t (B.2) 

where the linear operator L is given by 

£ ~5 ('-"5) < B ' 3 > 

The operator L acting on functions g that vanish in and <j) m is self-adjoint (with 
respect to the usual inner product (/,<?) = J^ m f(4>)g(4>) dip) and positive definite (to 
be precise L is acting on iJg (0, </> m )). In fact the Green's function associated to the 
operator L (i.e., the kernel of is explicitly given by e A ^ G(<p, <f>') where G is given 
by l|17|) . Notice that e A ^ G(</>, </>') is symmetric under the exchange of <f> and (f>' : and 
pointwisc positive in the interval (0, <p m ). 

Since the right side of (|B.2J| is nonegative for the solutions to (J2J) that we 
are considering, i.e., for the principal branch, and since L is positive definite, the 
corresponding solution g to the boundary value problem IjB.lfl is positive in (0, <p m )- 
By standard regularity properties of differential equations, the solution g is continuous 
in (0, 4> m )- In the neighborhood of the end points, the solution g of (|rj.l|) behaves as 
follows: g(<p) vanishes like \f^> near 0, whereas g((f>) vanishes like V </> m — <fi near <p m . 
To determine this behavior, notice that the solutions of (J2J corresponding to the 
principal branch satisfying the boundary conditions ©, Q are real analytic. Thus, 
in the neighborhood of x — (because of the boundary conditions), 4>{x) ~ <fi m — ex 2 , 
for some c > 0, hence, p — ~<p'(x) ~ 2cx near x — 0, and inverting (i.e., expressing 
x in terms of tfi we finally have p(4>) ~ y (fim — </>) in the neighborhood of <j> — <f> m . 
Analogously, p(4>) ~ yQ) in the neighborhood of <f> = 0. From the properties of the 
Green's funtion for the operator L (notice that G vanishes linearly near each of the 
endpoints) we see that the solution of l|B.l|l also behave as \ft$> near <fi = and as 
\J Wm ~ 0) near 4> = (f> m . 

Concerning the existence of the different integrals that appear on the right side 
of l|10(l when g is the solution of l|B.l(l . notice that h is continuous in (0,4> m ), and 
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although h is singular near the end points, h 1 / 3 behaves like 1/\Z4> (respectively like 
1 / V<Am — <t>) m the neighborhood of (respectively in the neighborhood of (j) m ), so 
ft, 1 / 3 is locally integrable near the end points. The other two integrals also exist, 
because g is continuous in [0,(/> m ]. 



Appendix C. 



In this appendix we derive the Euler-Lagrange equation for (|13|) . Consider the 
maximization of the functional 

-i 3 



K[g] 



h 1/3 ((j)) d<f> 



D 



where 



subject to 



W)=-~s" + «7+~s'. 



(C.l) 



(C2) 



g(4>)d<t> 



1. 



Introducing the constraint with a Lagrange multiplier A, wc obtain the Euler-Lagrange 
equation 

* h^dcf] 



1 



A d J_ l£_ 1 

A 1.2/3 1 _ O J».2 Vl 



/i 2 /3 2dy K h 2 / 3 ' 2dy 2K h 2 / 3 ' 



At the extrema h = 1/p 3 and the above equation becomes 

v 2 

A dp 2 1 d 2 p 2 ' 



P 



P - 77 



2 dy 2 dy 2 



+ A = 0. 



-Bct>+\ = 0.(C.3) 



(C.4) 



Identifying the Lagrange multiplier A with the eigenvalue r and observing that the 
integral in the equation above is equal to 1, we verify that the Euler Lagrange equation 
for extrema of K [g] is the original equation Q whose solution we seek. 



